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Abstract 

We argue that the entanglement entropy offers us a useful coarse-grained entropy 
in time-dependent AdS/CFT. We show that the total von- Neumann entropy remains 
vanishing even when a black hole is created in a gravity dual, being consistent with 
the fact that its corresponding CFT is described by a time-dependent pure state. 
We analytically calculate the time evolution of entanglement entropy for a free 
Dirac fermion on a circle following a quantum quench. This is interpreted as a toy 
holographic dual of black hole creations and annihilations. It is manifestly free from 
the black hole information problem. 
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1 Introduction 



The concept of holography, especially the AdS/CFT correspondence, has now been the 
most powerful tool to analyze quantum aspects of gravity [THS]. Since there have been 
no satisfactory and non-perturbative formulations of string theory, the direct calculations 
of quantum gravity effects often turn out to be quite difficult. However, the AdS/CFT 
equivalently transforms the problems into those in non-gravitational theories, typically 
large N gauge theories [3j. Since the latter ones are non-perturbatively well defined, we 
can in principle understand all quantum gravity effects at least in AdS spaces. 

In spite of enormous related developments, the AdS/CFT in time-dependent back- 
grounds have not been well understood at present. One of purposes of this paper is to 
understand better somewhat puzzling aspects of entropy in time-dependent backgrounds, 
which is also closely related to the black hole information problem pEJ. If we suddenly 
excite a CFT, it will eventually thermalize after a certain time. In AdS/CFT, this is 
dual to a black hole formation in AdS spaces [5l[6] . In a probe D-brane approximation, a 
thermalization is dual to a black hole formation in the induced metric on the brane [8]. 

In such a process, it is clear that the density matrix of the system is described by 
a time- dependent pure state in the CFT side, assuming that we start with a pure state 
before the thermalization. Therefore, the von-Neumann entropy is always vanishing. 
On the other hand, in the AdS gravity, since the black hole will be formed and the 
spacetime at late time is well approximated by the AdS Schwarzschild black hole, the 
entropy seems non-vanishing. For example, it is clear that the entropy which appears in 
the fluid description gets non- vanishing j9]. This might suggest an evolution of a pure 
state into a mixed state, which contradicts with the principle of quantum mechanics. 
To resolve this issue, we need to distinguish the microscopic von-Neumann entropy from 
a coarse-grained entropy. Only the latter can get non-vanishing in this thermalization 
process. We will argue that the entanglement entropy [T0HT2] offers us a convenient 
coarse-grained entropy in holographic setups. At the same time, we will show that the 
von-Neumann entropy of the total system is always vanishing in the black hole formation 
process employing the holographic calculation of the entanglement entropy [T3HT5] . 

From the CFT viewpoint, the black hole formation process is dual to a problem of non- 
equilibrium dynamics and therefore is quite complicated to analyze explicitly in general. 
Fortunately, recently there have been many progresses in understanding of properties of 
quantum field theories after thermalizations induced by instantaneous excitations called 
quantum quenches [T6H22] . We will consider a two dimensional CFT defined by a free 
Dirac fermion on a circle and calculate the time-evolution of the entanglement entropy 
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and correlation functions after a quantum quench, taking into account the finite size 
effect. Indeed, the behavior of the entanglement entropy we obtained implies that its 
gravity dual is supposed to be successive creations and annihilations of black holes in AdS 
spaces, though a free field theory in general corresponds to an extremely stringy region 
of the gravity side. This solvable example offers us a CFT description of a black hole 
formation and decay, which are manifestly free from the black hole information problem, 
because the system is clearly described by a pure state. Notice that in this problem, the 
standard thermal entropy is not useful as it is vanishing and the entanglement entropy is 
an alternative important quantity which can directly probe black holes. 

We would also like to mention that quite recently, remarkable numerical calculations 
of time evolutions of entanglement entropy have been done holographically in [23||24"]. 
When we essentially finished calculations for this paper, we found the paper [23], which 
has a partial but important overlap with ours about the argument which shows unitary 
evolutions of CFTs dual to the black hole formations by using the holographic entangle- 
ment entropy. A similar issue has also been discussed independently in the appendix of 
the paper [24]. 

The paper is organized as follows: In section two, we will explain a puzzle on the 
entropy in time-dependent backgrounds and resolve it by employing the entanglement 
entropy as a coarse-grained entropy. In section three, we will calculate the time evolving 
entanglement entropy in two dimensional CFT of a free Dirac fermion with finite size 
effect. We will interpret our results from the gravity dual viewpoint and relate them to 
the black hole creations and annihilations. In section four, we will study one and two point 
functions in our time-dependent CFT and confirm that they agree with our holographic 
interpretation. In section five, we will summarize our conclusion. In the appendix A, we 
present detailed calculations of two point functions in a free Dirac fermion theory on a 
cylinder using the boundary state formalism. 

2 Holographic Entanglement Entropy as Coarse-grained 
Entropy 

2.1 An Entropy Puzzle about Time-dependent Holography 

Consider holography in time-dependent backgrounds, especially the ones in AdS/CFT. 
We are particularly interested in a situation where a black hole horizon is created by a 
gravitational collapse. Such a time- dependent process has been considered to be dual to 
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thermalizations and indeed we can construct explicit examples in AdS/CFT setups [HE], 
where non-normalizable perturbations lead to a creation of hot plasma fluid in the dual 
Yang-Mills theory. Similar phenomena also occur in the probe D-brane setups regarding 
the thermalization as a formation of horizon in the induced metric [S] . 

An important characteristic physical quantity in time-dependent systems will be the 
entropy. Since the black hole is created after a certain time in our gravity background, it is 
quite natural that the entropy which is holographically calculated, gets non-vanishing at 
late time. However, in the dual CFT, the time evolution is described by a time-dependent 
external force in quantum many-body systems and thus the density matrix ptot(t) of the 
total system will evolve in a unitary way p{t) = U(t,t )p(t )U(t,t )~ 1 . Since the entropy 
of the total system is defined by the von-Neumann entropy 

S tot (t) = -Tr[p(t) logp(t)] = S tot (t ), (2.1) 

which is clearly time-independent. Therefore if we start with a pure state, then the 
entropy is always vanishing at any time, even if the black hole horizon appears in the 
dual geometry. This observation seems to lead to a contradiction between the gravity and 
CFT at first sight. 

Another confusing point is that there is no known unambiguous definition of entropy 
S t ot in the time-dependent gravity duals. This is partly because there are two different 
notions of horizons: event horizon and apparent horizon. Currently, there are several 
evidences that the latter will be more appropriate to define entropy [I3J[25]. However, 
still we need to specify the choice of time slices to calculate the apparent horizon and 
thus there are infinitely many different definitions. Remember that in static spacetimes, 
the time slice is uniquely chosen and there is no ambiguity. 

A closely related issue in the dual CFT side is that we can define a non-vanishing 
entropy even for a pure state if we perform a coarse-graining of the given total system. 
If the total system consists of a gas and a heat bath, then the entropy of gap is obtained 
by tracing out the heat bath and by the required coarse-graining due the fact that our 
observables are quite restricted compared with the microscopic degrees of freedom; see 
e.g. [26] for relevant quantum field theoretic formulations. Again this leads to infinitely 
many different definitions of (coarse-grained) entropy. Below we would like to resolve the 
above puzzle by paying attention to this coarse-graining procedure. 

The resolution of this puzzle is also closely related to that of the black hole information 
loss problem. This problem occurs since a massive object described by a pure state 
collapses into a black hole and finally ends up with a thermal gas generated by Hawking 
radiations, which looks like a mixed state. If this is true, a pure system has to evolve into 
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a mixed state after enough time and clearly contradicts with the principle of quantum 
mechanics. 

2.2 Entanglement Entropy as Coarse-grained Entropy 

Typically the coarse-graining is done by cutting out higher energy modes or higher order 
multi-particle interactions. However, this does not seem to allow straightforward calcula- 
tions in AdS/CFT setup. Therefore, here we would like to perform a coarse-graining by 
cutting out a spacial part of the total system. In this case, the obtained entropy coincides 
with the quantity called entanglement entropy. 

The entanglement entropy «SU(£) is defined for each subsystem A(t) of the total system. 
In quantum field theories, A(t) is specified by diving the total space at a fixed time t into 
two parts named as A(t) and B(t). Thus the number of choices of A are obviously infinite. 
The total Hilbert space H tot now becomes a direct product H to t = H A ® Hb, where the 
subsystem B is the complement of A. The reduced density matrix for A at any time t is 
given by 

p A {t) = Tr H B Ptot{t). (2.2) 
Finally we define the entanglement entropy as follows 

S A (t) = -Ti[p A (t) log p A (t)]. (2.3) 

This quantity is in general non- vanishing and a non-trivial function of time even for pure 
states. 

We can calculate Sa(^) via AdS/CFT by applying the holographic formula proposed 
in [Hj, which generalizes the holographic formula in static spacetime found in [13]. It is 
simply given by 

where Gn is the Newton constant of the AdS space. The surface 7a(£) is a codimension 
two surface in the AdS space which is defined by the extremal surface whose boundary 
coincides with the boundary of A(t). We require that 7a(£) is homotopic to A(t). If there 
are more than one such extremal surfaces, we pick up the one with the lowest area, as 
is required by the strong subadditivity of entanglement entropy [27]. Note that here we 
need to directly deal with Lorentzian spacetime without its Euclidean counterpart because 
it is time-dependent. Indeed, in [13], it was confirmed that such extremal surfaces are 
well-defined in important examples. Recently, this formula is applied to the AdS 3 Vaidya 
background and the conformal field theory results in [16] have been remarkably reproduced 
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in the work [23J. Moreover, this calculation has been extended to the higher dimensional 
cases in [24J. 

For our purpose, it is convenient to study a d + 1 dimensional time-dependent back- 
ground which is asymptotically global AdS^+i space, which includes a finite size effect 
because its boundary is R x S' d_1 . We assume that a black hole is produced at a certain 
time via a gravitational collapse and this process lasts only for a finite time. Our argu- 
ment below is heuristic and does not depend on the detailed form of the time-dependent 
solution. 

We calculate Sa at a time t when the black hole formation has been ended and when 
the spacetime is well approximated by a static AdS black hole. First we assume the 
length size of A (denoted by \A\ below) is smaller than the inverse temperature (3bh of 
the AdS black hole temperature. In this case, the extremal surface 74 (i) is localized near 
the boundary and surrounds the region A(t) [32]. Sa essentially consists of the area law 
divergent piece 

Area(&4) 

S * ^5=5— + ■-' ( 2 - 5 ) 

which are independent of the temperature. Notice that the condition \A\ << \B\ assumed 
here physically means that the system is heavily coarse-grained. As \A\ exceeds (3bh, Sa 
gains a finite and extensive contribution, which is essentially the thermal entropy for 
the subsystem A, in addition to the divergent part Sdi V given by (12. 5p . This thermal 
contribution comes from the part of the surface 7^ (t) which is wrapped on a part of the 
apparent horizon of the AdS black hole and thus is proportional to \A\ as in FigJH see [TITH^ 
[23j[2l] for explicit confirmations. An important fact is that though the apparent horizon 
requires the choice of the time slice, our extremal surface ta(0 is uniquely determined by 
the boundary condition. 

When this reaches the point \A\ = \B\, the situation begins to change importantly. 
If we naively speculate the surface ja (t) in a continuous way until \A\ exceed \B\, we 
may think that it will wrap on a more than half of the apparent horizon as in the right- 
up figure of FigJTJ However, since the horizon disappears at an earlier time, we can 
smoothly deform this surface into the one which wraps the opposite part of the horizon. 
Therefore the result of Sa as a function of \A\ becomes symmetric with respect to the 
point \A\ = \B\ as in the right-down figure of FigJTJ In other words, we can conclude that 
in this time-dependent background the entanglement entropy satisfies 

Sa = S B , (2.6) 

which is indeed known to be true when the total system is a pure state and thus agrees 
with our holographic setup. This sudden decreasing of the entanglement entropy occurs 
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because the subsystem B which we are tracing out gets smaller than the half of the total 
system and the information lost by the coarse-graining starts recovering^]. 

In the end, we can holographically calculate the von-Neumann entropy for the total 
system as follows 

S tot = lim (S A - S B ) = 0. (2.7) 

Thus, in the gravity side, we confirmed that the total entropy is always vanishing in spite 
of the black hole formation. As is clear from the previous argument, the essential point 
is that the horizon is time-dependent and gets vanishing at earlier time and the minimal 
area principle (among extremal surfaces) prefers the minimum i.e. Stot = at any time. 
In summary, in the time-dependent background which describes a thermalization of a 
pure state, the von-Neumann entropy of the total system is indeed vanishing also in the 
gravity side. This contrasts strikingly with the setup of an eternal AdS black hole [2H] as 
it is dual to a mixed state in a thermal CFT and has a non-vanishing thermal entropy, 
which precisely coincides with the value lim[B[_>.o (Sa — Sb) [30J . 

The non-zero entropy is obtained only after a certain coarse-graining. Among various 
such entropies defined in the dual CFT, only the entanglement entropy has its clear 
holographic dual in time-dependent backgrounds at present. From this viewpoint, we 
can define a coarse-grained effective entropy S e ff which looks analogous to the thermal 
entropy by 

S eff = 2(s A -Sf) . (2.8) 

V / \A\ = \B\ 

We subtracted the entanglement entropy before the thermalization. This clearly 
includes the divergent part given by (12.51) and therefore S e ff is finite. The reason 
why we put the factor two is A covers only a half of the total system. We would like to 
argue that (12. 8 j) is a definition of coarse-grained entropy which is uniquely calculable both 
in gravity and CFT side of the AdS/CFT. Notice that though the calculation of entropy 
from an apparent horizon requires the choice of time slice and is ambiguous, our entropy 
S e ff has no ambiguity as the extremal surface condition of ^y A (t) determined where the 
surface wraps the apparent horizon. 

Finally, we would like to come back to the relation to the black hole information 
problem pE]. From the above argument, it is clear that there is no information loss in 
the gravitational collapsing process because S to t = is always satisfied. However, we 
cannot study the most important process of the decay due to Hawking radiations because 



3 A similar behavior has been known in the analysis of quantum information in evaporating black 
holes [28| . This has been considered to be crucial to resolve the information paradox of quantum black 
holes because the information may be recovered after more than half of a black hole has been evaporated. 
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Continuous Deformation 

Black Hole Formation in AdS 




Figure 1: The minimal surface 7^ for the holographic calculation of entanglement entropy 
in the AdS black formation (left). Though for an eternal AdS black hole we have ia 7^ 1b 
(upper right), in our case of black hole formation, we actually find that ^a — 7s (lower 
right). This is because the horizon vanishes at early time. 

the absence of the holographic formula of entanglement entropy beyond the supergravity 
approximation. Therefore, for this purpose, it is a better idea instead to study its CFT 
dual, directly. We will analytically investigate such an example in the next section. 

3 Time-dependent Entanglement Entropy in 2D CFT 
with finite size via Quantum Quench 

In the previous section, we studied the behavior of coarse-grained entropy by using the 
holographic entanglement entropy. In this section, we would like to complete the argument 
by presenting an explicit calculations with a finite size effect in two dimensional CFT and 
by showing that the results agree with the holographic ones. We would like to analyze 
how the entanglement entropy evolves in time-dependent backgrounds which describe the 
thermalization phenomena. At the same time, this offers us a toy gravity dual of a black 
hole formation and evaporation, where there is no information loss problem and where 
we can solve the model exactly. 

We will particularly look at a specific class of time-dependent system which is called 
quantum quench, whose field theory descriptions have been given in [T6H22] . The quantum 
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quench is a phenomenon when we suddenly change parameters of a given system such as 
masses and coupling constants. The evolutions of entanglement entropy under quantum 
quenches have been studied in [l6j[T9j[20] for infinitely extended systems. 

To make an exact calculation possible, we consider a free Dirac fermion theory in 
two dimension. In free field theories, the thermalization via the quantum quench is not 
a standard one described by the (grand-) canonical ensemble because the different mo- 
mentum modes are decoupled and do not mix [T8l[2TJ . Notice also that the free CFT 
corresponds to string theory with a lot of quantum corrections via the AdS/CFT. Nev- 
ertheless, this system is enough for our purpose qualitatively as we will see later. As far 
as the authors know, our result will be the first analytical calculation of time-dependent 
entanglement entropy with a finite size effect. The calculation of entanglement entropy at 
finite temperature and with finite size has been calculated in [30] for a free Dirac fermion. 

3.1 How to Calculate Entanglement Entropy under Quantum 
Quench 

The setup of quantum quench that we have in mind is as follows. Start with a Dirac 
fermion theory with a mass m in two dimension. Then at the time t — 0, suddenly we 
turn off the fermion mass. Then the system after this quantum quench i.e. t > will be 
described by a time-dependent background in the c = 1 2D CFT defined by a massless 
Dirac fermion whose Hamiltonian is written as H. If we denote the state just after the 
quench by l^o); then the density matrix at a time t > is given by 

Ptot(t) = e- itH \y )(y \e UH . (3.9) 

To find |\l/o) we would like to follow the prescription in [16J. We expect that in the infrared 
limit |\l/o) should flow into a boundary fixed point. Therefore we can approximate |\l/o) 
by a boundary state \B) with the UV modes filtered out 

\V Q )~e- eH \B) } (3.10) 

where e is the parameter of UV filter and is inversely proportional to the original fermion 
mass e ~ m~ l . If we denote the UV cut off or lattice spacing of this field theory by auv-, 
we need to take e >> auv, which is clearly satisfied in the continuum limit auv — > 0. 
In this way we find that the reduced density matrix is eventually written as 

Ptot(t) = e- UH -' H \B}(B\e UH ^ H . (3.11) 
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In the path-integral description [IB] , the spacetime where the Dirac fermion lives is 
the cylinder geometry given by the direct product of an interval (i.e. Euclidean time) and 
a circle (i.e. space). The length of the interval is 2e and the one of the circle is normalized 
into 2n. We introduce a complex coordinate (y, y) for this spacetime such that 

y = T -ia, y = r + ia, (0 < r < 2e, < a < 2vr). (3.12) 

The Dirac fermion on this spacetime is decomposed into the left-moving part 

{ip L {y),i>L{y)) and the right-moving part (ip > R (y) ,ip > R (y)) . 

We define the subsystem A(t) as an interval in the limy direction at time e + it. We 
can parameterize the location of the two ends points of this interval as 

(yi,y~i) = (e + it + iai,e + it -iai), and (2/2,2/2) = (e + it + i(T 2 , e + it - ia 2 ). (3.13) 

In the path- integral description [IH1ES], the reduced density matrix pA(t) is obtained 
from ptotif) (13.111) by gluing the subsystem B and the two boundaries on the subsystem A 
are left as they are (see the left figure in Figj2]). Then the important quantity Trf/j^t)^] 
can be found as the partition function of the free Dirac fermion on a iV-sheeted Riemann 
surface S^r [35] as shown in the right figure in FigfJl is defined by the N copies of 
the cylinders (13. 12ft which are glued successively along the subsystem A. This quantity 
is used to define so called the Renyi entropy. Recently, remarkable properties of Renyi 
entropy and their consistency with holographic results have recently been studied in [31] . 

In this way, the entanglement entropy can be computed via 

S A (t) = -^\ogTr[ PA (t) N }\ N=1 . (3.14) 

By employing the replica trick, we can regard the system is given by N free Dirac 
fermions (^ a \^) a = 0, 1, • ■ ■, N- 1 G Zn on a single cylinder with the twisted bound- 
ary conditions at the two points ip^ \e 2m y\) = il)^ +1 \yi) and ^\e 2m y2) = ip^ 1 \y2)- 
Since the lagrangian is gaussian, its form remains unchanged even after we perform the 
discrete Fourier transformation vjj^ — > ^= YltLo 1 e^^tp^ [32] . The interactions be- 
tween N fermions only appear through the twisted boundary condition which are now 
diagonalized as follows 

= e^\ Vl ), 4V^ 2 ) = e-^f{y 2 ). (3.15) 

These boundary conditions can be equally replaced with the twisted sector vertex opera- 
tors with the lowest conformal dimensions. Finally, we can calculate the entanglement 



9 




n-sheets 



Ze 




Tr[(p A ) n ] 



Figure 2: The path-integral calculation of the reduced density matrix [pa] a/3 (left) and 
the trace Trp A (right). 



entropy from the formula 

Tr\p A (t) N ]= J] (^ a \y 1 ,y 1 )^- a \y 2 ,y 2 )) cylmder 

(B\e- 2£H a^(y 1 ,y l )a^(y 2 ,y 2 )\B) 



N-1 

2 



N-1 
2 

n 



(3.16) 



JV-l 
2 



(B\e~^ H \B) 

Here we normalized the two point functions such that Tr[p^(t)] = 1. 

3.2 Bosonization and Two Point Functions on Cylinder 

A useful fact of the free massless Dirac fermion (ip, ip) is that it can be bosonized into a 
free scalar X via the standard formula 

My) = e lXL{y \ My) = e- lXLiy \ My) = e" [ " m , My) = e~ lX ^\ (3.17) 

Our normalization of the two dimensional CFT in this paper is a' = 2 in string theory 
world-sheet |33j and the free Dirac fermion is equivalent to the compact boson X at 
the radius R = 1. If we write the primary fields as V(k L ,k R ){y, y) — e %kLXL ^ +lkRXR ^\ 
then there are two candidates of the twisted vertex operators, which lead to the twisted 
boundary conditions (I3.15p . 



<?i\y,y) = V(« _ »)(y,y) 



ijj{X L {y)+X R (y)) 



(3.18) 
(3.19) 
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which have the lowest dimensions if a runs a = — — ^j^, • • •, ^jr- [30] . 

What we need to calculate to obtain the entanglement entropy are the two point func- 
tions of ( 13.181) or (13.191) on the cylinder. This can be found by performing calculations 
using the explicit form the boundary state \B) (see e.g. |43j) as summarized in the ap- 
pendix A. One may think that there are two choices i.e. the Neumann and Dirichlet 
boundary condition in the free scalar theory. If we returns to the Dirac fermion theory, 
the Neunmann boundary condition relates ipL and ipL to ipR and ipR at the boundary, 
while the Dirichlet one does i/jl,r to ipR,L, respectively. We can see that the discrete 
Fourier transformation, employed to map the replicated fermions into the new ones which 
satisfy (13.151) . is consistent only when we choose a[ a ^ or a 2 for the Neumann or Dirichlet 
boundary condition, respectively. 

The final result of the normalized two point functions in both Neumann and Dirichlet 
is given by the same expression 

a 2 

/2ieN6 \n I t+it _|_ _£_\2ie\nn I e+it a_ I 2ie\ I \ 

1^(^1216)12.^(^1^)12 J ' ^> 

where we defined a = a 2 — <J\ ■ Refer to appendix A for a derivation of this result (I3.20|) . 
About the definition of the eta function r/(r) and theta function 9i(v\t) we followed the 
convention in 1331. 



3.3 Time Evolution of Entanglement Entropy with Finite Size 

By substituting the two point functions (13. 20f) into (13. 16ft and (13.141) . we can calculate the 
entanglement entropy ^(t). Using the formula 

iV-l 

a 2 1 I 1\ , 

£ W = u\?-Ti)> (3 ' 21) 

0= — 

we obtain the entanglement entropy 

1 \Q ( JL 1 2il ) 1 2 \Q I gjeN 12 

S A (*, a) = - log — - ln \^ VL ^^\ l2l ^ _ 2 , (3.22) 



where we recovered the cut off aw dependence in an obvious way. Notice that this is 
independent from the choice of the boundary condition (Neumann or Dirichlet one). From 
this it is clear that SU(£) satisfies 

S A {t, a) = S A {t, 2n-a) = S B (t, a). (3.23) 
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This is consistent with the assumption that the total system is a pure state and the 
von-Neumann entropy for the total system is vanishing. 

After the modular transformation, we can rewrite ( 13.221) as follow£| 

q (+ x _ 1 1 2e 1 , l^i(flg)l 2 -l^i(^rli)l 2 ,„ 0A , 

The explicit form of S A (t, °") is plotted in Figj3]as a function of a and FigJUas a function 
of t. In the limit e — >■ 00, where the quench disappears, we reproduce the standard 
result [SUES] 

S"a = - log f — sin - j . (3.25) 



3 V ^^ 2 

In the opposite limit e — > 0, which corresponds to the infinitely extended space, we can 
easily reproduce the known result in [16] at the central charge c = 1 

S A = S div + ^- (0<t<^) 
be 2 

= S dm + ^ (*>f), (3.26) 

where we separate the divergent part S^w = f log ■ The physical meaning of the 
constant term can be qualitatively understood as that of a thermal CFT gas with 
length a and at the effective temperature 

Teff = \ e (3.27) 

We can also calculate the effective coarse-grained entropy introduced in (12. 8p as follows 

9 f0l— ) \9 (siiil^ii \\ 

s eff (t) = 2 (s A (t, rc) - s A (o, tt)) = - log ^Tk^VjAri^ ' ( 3 - 28 ) 

which is clearly finite and does not depend on the cut off auv- Then we notice a periodicity 
in time 

S A (t, a) = S A (t + vr, a), S eff (t) = S eff (t + tt), (3.29) 

as is clear in in FiglU This short periodicity is peculiar to our free CFT, where excitations 
triggered by the quench propagate at the speed of light and return to exactly the same 
state after it goes half of the circle. This is because the energy of oscillators included in 
the boundary state is always an even integer. In generic field theories this does not occur. 
However, when the size of the space manifold is finite, then the Poincare recurrence leads 



4 Actually, in this expression, the values of Ox | § ) , Ox ( ^ + f§ | §) and Ox ( ^ - g | are always 
all real even before we take their absolute values. This justify our analytical continuation Re[y] e + it. 
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Figure 3: The plot of 5^(^/2, a) — 5^(^/2, 0) as a function of a at e = 0.2. 

to a much longer periodicity of time which is estimated as an exponential of a coarse- 
grained entropy At ~ e Seff , which has been the crucial fact to resolve the information 
problem for eternal AdS black holes [29 |[37j H0]. 

The AdS/CFT implies that a holographic dual of gravity on AdS3 with significant 
stringy corrections is a weakly (or free) coupled field theory^, which shows the short time 
periodicity like (13.291) . In the gravity side, this oscillating behavior of S^t) and S e ff(t) 
can be interpreted as successive creations and annihilations of black holes as sketched in 
Figj5j In this way, our free field theory example offers us a holographic dual of black hole 
formations and evaporations in an extremely stringy region and this is manifestly free 
from no information loss problem as the total system is described by a pure state. Even 
though our argument on the information problem is along the line with [28 | |29 |[371 - H0] . 
our approach is new in that we explicitly treat a real time-dependent process of black 
hole formations by employing its holographic dual, rather than focusing on the static 
black holes. Notice again that in this example, the size of black holes are probed by the 
entanglement entropy and the standard total von-Neumann entropy is vanishing and not 
useful. 



4 Time-dependent Correlation Functions 

In previous section, we observed that S e ff is oscillating in our time-dependent system. We 
then identified the dynamics in CFT is dual to successive productions and annihilations 

5 Strictly speaking, in the standard examples of AdS3/CFT2, we also need to take into account the 
free scalar field sector in addition to the fermion sector we analyzed. However, since both are free field 
theory, we expected that the qualitative result will not change. 
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Figure 4: The plot of S eff {t) = 2 (S A {t, n) - S A (0, vr)) as a function of t at e = 0.2. 



Time 




Quench 



Figure 5: Quantum black hole creations and annihilations in AdS space by the quantum 
quench as obtained from FigJU 
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of quantum black holes. We can naturally imagine the mechanism of its dynamics as 
follows. First, a black hole is created by a quantum quench and later it evaporates via a 
quantum analogue of Hawking radiations, which are usually semiclassically analyzed. But 
the radiations are completely reflected by the time-like boundary of the AdS. Therefore 
they again concentrate at the center of AdS and recreate a black hole. To confirm this 
idea, we would like to find time-dependent signals of radiations]^ Though we will perform 
explicit computations for a free Dirac fermion below, we can equally apply the results in 
this section to correlation functions in a free scalar theory compactified at R — 1. 



4.1 One Point Functions 

Assume that an operator O in CFT is dual to a field ip in AdS. The bulk to boundary 
relation in AdS/CFT [UJ argues that the one point function (0(x)) in the CFT is pro- 
portional to the value of the field tp near the AdS boundary. Since we are not deforming 
the CFT itself, the profile of <p is supposed to be normalizable. 

The most important one point function in time-dependent backgrounds will be that 
of the energy stress tensor. Due to its conservation law, this clearly becomes time- 
independent. It is straightforward to calculate the energy density £ n of the left-moving 
oscillator a n from (13 .lip as follows 

1 Ti 

£ » = to {a - nan) = 2«(<*»-iy (430) 

This indeed obeys the Bose-Einstein distribution with the effective temperature estimated 
in f)3.27p . This temperature can also be seen in the reduced density matrix when we trace 
out the right-moving oscillators a n (or equally left-moving ones) 



PL = Tr R e-* H \B)(B\e-* H = FT Ve" ta ^=^| ><0|^L, (4.31) 

■ni — 1 1 <r> — ( 1 v v 



m=0 n=0 



where we omit the zero modes. The total energy density in the e — > limit can be found 

as usual £ tot ^ 2 J2n=i - f T eff 

Even though we confirmed that our system is thermalized, we cannot find signals of 
black hole creation and annihilation from the time- independent one point function f)4.30p . 
Motivated by this we would like to calculate another one point function which is not dual 



6 For asymptotically flat black holes, we can successfully calculate radiations from thermalized D- 
branes by taking into account the coupling between open string fields and close string fields as in |36j . 
Below we would like to consider time-dependent radiations from black holes in AdS spacetime rather 
than flat spacetime. 
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Figure 6: The plot of the one point function e tkX ^' ^ as a function of t at e = 0.2 and 
fc = 3. 



to any conserved quantities. Especially we consider the one point function ( e lfcX (*>°")) for 
the Dirichlet boundary state for the initial state l^o); where k can take only integer values 
as we set R — 1. After a calculation similar to the appendix A, we finally obtain (we set 
a = 0) 

3 



JkX(t t Q) 



k{e+it) | 2ei 







3k 2 


101 


fe+it 

V 7T7 


7T / 


fe 2 



(4.32) 



As is clear from the plot in FigJHl it has peaks at t — 0, 7T, 27r, ■ • - . This behavior is actually 
natural from our interpretation in terms of black hole creations and annihilations. As the 
result of the entanglement entropy suggests, a black hole is created at t — and reaches 
its maximal size at t = |. Therefore, we expect that the radiations from the black hole 
will be the most strong at t = |, • • •, remembering the periodicity n. They will reach 
at the AdS boundary after the propag 

ationtim.fi A t = V2. This e^ains the peaks of 
(e jfcX (*'°)) since the one point function is dual to the value of bulk scalar field near the AdS 
boundary and its square should be proportional to the strength of radiations. Notice that 
this non-vanishing one point function is peculiar to time-dependent black holes. This is 
because the one point function, which is holographically dual to a classical radiation, is 
vanishing in the static thermal CFT. 



7 This can be found from the metric of the global AdS3 : ds 2 — — cosh 2 pdt 2 + dp 2 + sinh 2 pdO 2 as 
At = co ^ = § • In the presence of black holes one may think that the propagation time of a massless 
particle will be changed. Indeed, in a semiclassical black hole, this happens and affects the structures 
of singularities of two point functions |42) . However, we do not consider this modification as we are 
interested in a region in gravity where quantum corrections are so large that the even horizon does not 
seem to be defined using a metric. 
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4.2 Two Point Functions 



The another reason that we are interested in time dependent correlation functions is 
because we would like to discuss the information loss problem of black holes [4J. In 
[211133110], it has been argued that information loss problem is dual to the quasinormal 
behavior (exponential decay) of time dependent two point functions in a large N gauge 
theory at finite temperature: 



(0{h)0{t 2 )) ~ expf-cjt! - * 2 



(4.33) 



where c is a positive real constant. These behavior breaks the unitarity of the theory 
because this shows that any fluctuations that we add to the initial state will eventually 
vanish. In the same papers, it has also been conjectured in a finite N gauge theory, the 
decay of the correlation stops at the value of order exp[— 0(N 2 )], then correlation grows, 
and repeats this behavior because of the Poincare recurrence of the theory. 

Indeed, this quasinormal behavior can be seen in our model. If we set y\ = e + it, 
1)2 = e and ki = kn = k in the two point functions ( 1A.44I) for the Neumann boundary 
condition, we find 



G(t) = (expikX(t,0)exp-ikX(0,0)) 



N 



2ie 
V 1 — 

71 



OA - 



01 



(g+g) | 2ie 



a I t I 2ieV 
01 \2^\—> 



n ( 2e+it I 2ie \ ■ 
^l \-2ri~\ IT) 



A- 2 



(4.34) 



In FigJTJ we show the plot of G(t). It has the periodicity 2tt as a signal propagates 
from cr = and it has to go back to the same point. The divergences at t G 2nZ come 
from the usual short distance behavior of the operator product. If e t <C 1, we can find 
the following exponential behavior of G(t) 



G(t)=;/(e)exp 



2+1 



Tlk t 

'~2T 



(4.35) 



for a certain time-independent function /(e). Due to the 2n periodicity, after G(t) gets 
decreased exponentially, it reaches its (non-zero) minimum at t — n and then again begins 
to increase. 

Finally, it is also intriguing to study the time evolution of the two point function with 
separated points (expikX(t, a) exp — ikX(t, 0)) for Neumann boundary condition. If we 
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Figure 7: The plot of G(t) as a function of t at e = 0.2. 
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Figure 8: The plot of F(t, a) as a function of t at e = 0.2. 
use the formula flA.441) of Appendix A, and set y\ = e + it + zer and ?/2 = e + it, we have 



F(t,cx) = (expifcX(t,<r)exp-ifcX(t,0)) 



2ieY 



iV 



V — 



I Q ( I 2ie \ 



\2k 2 



71 J 



\n ( a I 2ie\ 



\2k 2 



2(e+it)+ia i 2re 

27T2 ' 7T 



fc 2 



2(e+^)-ia- 1 2it 
2ni I 7T 



A- 2 ' 



In FigJBJ we plot F(t, a). The graph shows that the correlation becomes strong when 
S e ff is large (compare with FigH]). This can be regarded as a further evidence of the 
successive productions and annihilations of quantum black holes in our system. 
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5 Conclusions 



In the first part of this paper, we pointed out a puzzle. Because a thermalization process 
in the CFT side should be unitary, the von-Neumann entropy of the final state must 
remain vanishing. On the other hand, in its gravity dual, we seem to have non-vanishing 
entropy as the final state includes a black hole. These observations look contradict with 
each other. However this is not the case. We showed that even in the gravity side, the von- 
Neumann entropy of the final state is zero by using the property Sa = Sb of holographic 
entanglement entropy for an arbitrary pure state. We proved this by making use of the 
fact that the apparent horizon of the black hole vanishes in the past. We then proposed 
a direct generalization of thermal entropy by using the entanglement entropy. From 
a holographic point of view, this quantity seems reasonable because S e ff is essentially 
proportional to the area of the apparent horizon of the corresponding black hole. 

Nevertheless, if quantum corrections are taken into account, one might think that the 
entropy of the final state in the gravity dual can be also nonzero due to the Hawking 
radiation. This observation is closely related to the information loss problem. In the 
latter part of this paper, we calculated the time evolution of entanglement entropy fol- 
lowing a quantum quench. Especially, we considered two dimensional Dirac fermions on 
a cylinder. The result shows that the von-Neumann entropy is always vanishing and the 
coarse-grained entropy S e ff is oscillating as a function of time. This offers us a clear toy 
holographic dual of a formation and decay of black holes, which is free from the informa- 
tion problem. We interpreted this fact as the successive productions and annihilations of 
quantum black holes. We also studied the behavior of one and two point functions and 
confirmed that they support this interpretation. We found that the two point functions 
cannot be zero and is periodic in time which shows some recurrence of our free system. 

Acknowledgments We are grateful to M. Headrick, W. Li, S. Ryu, Y. Srivastava and S. 
Sugimoto for useful comments. TT and TU are supported by World Premier International 
Research Center Initiative (WPI Initiative), MEXT, Japan. The work of TT is also 
supported in part by JSPS Grant-in-Aid for Scientific Research No. 20740132, and by 
JSPS Grant-in-Aid for Creative Scientific Research No. 19GS0219. TT is supported by 
World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. 



19 



A A Derivation of Two Point Functions on Cylinder 



Here we would like to present an explicit derivation of two point functions (I3.20p on 
cylinder. We work in the description of a massless scalar field X via the bosonization of 
the Dirac fermion. We set a' = 2 in string theory convention and follow [33] about the 
notations of the mode expansions of a scalar field and the definitions of theta functions. 

We employ the complex coordinate (y,y) on the cylinder defined by ( 13.121) . In the 
operator formalism, the scalar field is quantized as usual 



X L = x L - ip L y + % } — e 



® r " —my 



m 



X R = x R - ip R y + i V ^ e ~ m \ (A.36) 
L — ' m 



where the commutation relations are given by 

[xl,Pl] = h [ x R,Pn] = h [oc m ,a n ] = md n+mfi} [a m ,a n ] = m5 n+mfi . (A.37) 

At the compactification radius R, the momenta are quantized as follows 

n wR n wR 

Pl = 5 + — , P* = s -— , (A.38) 

in terms of the integers n (momentum) and w (winding). The (un- normalized) two point 
function of the normal ordered vertex operator Vi^ L ,h R ) ='■ e lkLXL+lkRXR : is written as 
follows 

( V (k L ,k R ) (Vl, yi)V(-k L ,-k R )(y2, i/2)) cylinder 

= (Sie- 2 ^^,^)^,^)^,^)^,^)^}. (A.39) 
The Hamiltonian is given in terms of the Virasoro generators as H = L$ + Lq — r=. 



A.l Neumann Case 

The boundary state \B) for the Neumann boundary condition [13] is given by 

oo 

\B) N =M e -^n=lW«-» 5 -n \W), (A.40) 

W= — OO 

where M is the normalization factor whose explicit value is not necessary for our purpose. 
For bosonization procedures in the boundary state formalism, refer to [H]. 
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The zero-mode part of flA.391) can be calculated by using the Baker-Campell-Hausdorff 
(BCH) formula as follows 



E 

w=— oo 
oo 

= E 



/ w l e ~2eH e ik L (x L -ip L yi)+ik R (x R -ip R y 1 ) -ik L (x L -ip L y2)-ik R (x R -ip R y2) 



w) 



R w g _R 

2 e 2 



(k L w(y 1 -y2)-k R w(y 1 -y 2 )) e ^r{yi-y2)+^ L (y2-yi) 



(A.4i; 



W) = — oo 



The massive modes can be computed by employing the BCH formula, leading to the 
identity 

1 



^ e -a/3z e a L a+a R /3 e b L a++b R P+ e ~a+ /3+ |g^ 



1 - Z 



a L b L+ a R b R- a L a R- zb L b R 



(A.42) 



where we assume the commutation relation [a, a + ] = \J3, = 1. Especially we can 
perform the calculations for massive parts by taking z = e~ 4ne and 



' e -nyi _ e ~ny2 s 



n 



n 



(e nVl -e ny2 ), 



)i a R ~ 
k R 

m = 



n 



in (1A.42|) . By expanding (1 — e 



4ne\ — 1 



n 



00 „-4mne 



(A.43) 



Eoo 
m=0 



and changing the order of the 



summation w.r.t n and m, eventually we can factorized the result into 9 functions. The 
final result is given by 

(B\e- 2eH V {kLM) (y 1 ,y 1 )V^ kL ^ kR) {y 2 ,y2)\B) N 



M 2 



e 2 e 2 

w=— OO 
2ie\ 3 



(k L w(y 1 -y 2 )-k R w(y 1 -y2)) 



2ie\ 3 
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V 27TI 
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» 27TJ 
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k ) V1 \ 2ni 


2ie\ 

K ) 



khk R 



(A.44) 



After we substituting the values ( I3.13P and the twisted vertex operators ( I3.18P at the free 
fermion radius R = 1, we find that the second line of the above expression is rewritten as 



TZ—e-v _e a (offl + th (o|f) 



vffl 



vffl 



(A.45) 



where the first and second term are interpreted as the NS and R-sector of the Dirac 
fermion. Therefore it is equal to (B\e~ 2eH \B) and can be neglected in order to find the 
normalized two point functions as we did in f )3.16p . In this way we obtain the result (13.201) . 
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A. 2 Dirichlet Case 

For the Dirichlet boundary condition, the boundary state is given by 



\B) D =W e E "=i 



-nw—n 



By repeating similar calculations, in the end, we find the two point functions 

(B\e~ 2eH V( kLikR )(yi, y 1 )V^k L ,-k R )(y2, V2)\B)d 



N 2 



E 



2 

g~ ~gr e f ( k L (yi - V2 ) +k R ( yi - y 2 ) ) 



, - 1 - V 27TJ 1 7T / / \ V 2-7TI 1 7T / 



0i (V 51 

x V 27TJ 
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TT J M 2m 


2ie\ 
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/i (yi+m 
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2ie\ 
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(A.46) 



-k L k R 



After we substituting the values (13.1 3D and the twisted vertex operators (13. 19ft at the free 
fermion radius R = 1, we find that the second line of the above expression is rewritten as 



yc 



-2n 2 e 



^ (oiv) 



(A.47) 



which includes only the NS-sector. 
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